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£N| ■ We present an extension of the Monte Carlo Event Generator WHIZARD which includes 

chargino production at the ILC at NLO. We include photons using both a fixed order 
■ and a resummation approach. In the latter, leading higher order corrections are auto- 

X^/^ ' matically included. We present results for cross sections and event generation for both 

methods [TJ.This is an updated version of the results presented in [2]. 



^ ■ 1 Introduction 
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In many GUT models, the masses of charginos tend to be near the lower edge of the supcr- 
'■ partner spectrum, and they can be pair-produced at a first-phase ILC with cm. energy of 

, 500 GeV. The precise measurement of their parameters (masses, mixings, and couplings) 

is a key for uncovering the fundamental properties of the MSSM [3J. Regarding the exper- 
, imcntal precision at the ILC, off-shell kinematics for the signal process, the reducible and 

' irreducible backgrounds [4], and NLO corrections need to be included. We here present the 

inclusion of NLO chargino production where corrections can be in the percent regime. 

(N 
00 ■ 

2 Chargino production at LO and NLO 

The total fixed-order NLO cross section is given by 

(7tot(s,TOe) = CTBorn(s) + (T v+S (s, A£ 7 ,m^) + <7 2 ^ 3 {s, A£ 7 ,m^), (1) 



where s is the cm energy, m e the electron mass, and A E 7 the soft photon energy cut 
dividing the photon phase space. The 'virtual' contribution a v is the interference of the 
one-loop corrections [5] with the Born term. The collincar and infrared singularities are 
regulated by m e and the photon mass A, respectively. The dependence on A is eliminated 
by adding the soft real photon contribution a s = / so ft fBorn(s) with a universal soft factor 
/soft(— r 1 ) 0. We break the 'hard' contribution <T2^3(s, A£? 7 , ml), i.e., the real-radiation 
process e~e + — ► X^Xj'li hito a collinear and a non-collinear part, separated at a photon 
acollinearity angle A# 7 relative to the incoming electron or positron. The collincar part is 

2 

approximated by convoluting the Born cross section with a structure function /(x; A# 7 , — =■) 
[7] . The non-collinear part is generated explicitely. 

The total fixed order cross section is implemented in the multi-purpose event generator 
O'Mega/WHIZARD [3 [9] using a 'user-defined' structure function and an effective matrix 
element |7W e fr| 2 which contains the Born part, the soft-photon factor and the Born-1 loop 
interference term. In the soft-photon region this approach runs into the problem of negative 
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Figure 1: 9 -dependence of effective squared matrix element (i/s = 1 TeV).Left figure: fixed 
order effective matrix element; right figure: effective matrix element with the one-photon ISR part 
subtracted. Solid line: Born term; dashed: including virtual and soft contributions for A_B 7 = 
10 GeV; dotted: same with A£ 7 = 0.5 GeV. A0 7 = 1°. 



event weights [TU]: for some values of 8, the 2^2 part of the NLO-corrected squared matrix 
element is positive definite by itself only if AE 7 is sufficiently large, cf Fig.[TJ To still obtain 
unweighted event samples, an ad-hoc approach is to simply drop events with negative events 
before proceeding further. 

Negative event weights can be avoided by resumming higher-order initial radiation using 
an exponentiated structure function /isr [H]. In order to avoid double-counting in the 
combination of the ISR-resummcd LO result with the additional NLO contributions |5], 
we have subtract from the effective squared matrix element the soft and virtual photonic 
contributions that have already been accounted for in er s+v . This defines lA^Jj^l 2 = |-M ff| 2 — 
2/ so ft,iSR | Bom | 2 which is positive for even low A£" 7 cuts for all values of 9 (cf Fig. QJ, such 
that unweighting of generated events and realistic simulation at NLO are now possible in all 
regions of phase-space. Convoluting this with the resummed ISR structure function for each 
incoming beam, we obtain a modified 2 — > 2 part of the total cross section which also includes 
soft and collincar photonic corrections to the Born/one- loop interference. This differs from 
the standard treatment in the literature (cf eg. [5]) where higher order photon contributions 
are combined with the Born term only ("Born+"). The complete result contains the hard 
non-collincar 2^3 part convoluted with the ISR structure function: 

d Xl dT 2 fjf R \x 1 )rif R \x 2 )\M^\ 2 + I d Xl dT 3 fjf R \x 1 )f^ R \x 2 )\M 2 - 3 \ 2 

JA(E,e) 

(2) 

3 Results 

Fig. [5] compares the A£" 7 dependence of the numerical results from the semianalytic fixed- 
order calculation with the Monte-Carlo integration in the fixed-order and in the resummation 
schemes. The fixed-order Monte-Carlo result agrees with the semianalytic result as long as 
the cutoff is greater than a few GeV but departs from it for smaller cutoff values because 
here, in some parts of phase space, |-M ff| 2 < is set to zero. The semianalytic fixed-order 
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Figure 2: Total cross section dependence on AE-y\ 'sa' (dotted) = fixed-order semianalytic result; 
'fix' (dashed) = fixed-order Monte-Carlo result; 'res' (long-dashed) = ISR-resummed Monte-Carlo 
result; (dash-dotted) = same but resummation applied only to the 2 — > 2 part. A9- 1 = 1°. LO: 
Born cross section. 



result is not exactly cutoff-independent, but exhibits a slight rise of the calculated cross 
section with increasing cutoff (breakdown of the soft approximation). For AE 1 = 1 GeV 
(10 GeV) the shift is about 2permil (5permil) of the total cross section. The fully resummed 
result shows an increase of about 5 permil of the total cross section with respect to the fixed- 
order result which stays roughly constant until AE 1 > 10 GeV. This is due to higher-order 
photon radiation. 

In Fig. [3] we show the binned distribution of the chargino production angle obtained using 
a sample of unweighted events. It demonstrates that NLO corrections (which, for total cross 
sections, are in the percent regime and can reach 20% at the threshold) are important and 
cannot be accounted for by a constant K factor. Figure 3] shows the magnitude of second 

and higher order photonic effects in different schemes. Resummation effects are clearly in 
the percent regime and cannot be neglected. For y/s > 500GeV, the convolution of the 
interference term with /isr additionally changes the sign of the higher order corrections. 
For more details, cf. [121 113j . 

4 Conclusions 

We have implemented NLO corrections into the event generator WHIZARD for chargino pair- 
production at the ILC with several approaches for the inclusion of photon radiation. A 
careful analysis of the dependence on the cuts A£ 7 , A 9 reveals uncertainties related to 
higher-order radiation and breakdown of the soft or collinear approximations. To carefully 
choose the resummation method and cutoffs will be critical for a truly precise analysis of real 
ILC data. The version of the program resumming photons allows to get rid of negative event 
weights, accounts for all yet known higher-order effects, allows for cutoffs small enough that 
soft- and collinear-approximation artefacts are negligible, and explicitly generates photons 
where they can be resolved experimentally. Corrections for the decays of charginos and 
non-factorizing corrections are in the line of future work. 
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Figure 3: Polar scattering angle distribution for an integrated luminosity of 1 ab _1 at *Js = 1 TeV. 
Left: total number of events per bin; right: difference w.r.t. the Born distribution. LO (dotted) 
= Born cross section without ISR; fix (dashed) = fixed-order approach; res (full) = resummation 
approach. Cutoffs: AE-y — 3 GeV and A# 7 = 1°. 
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Figure 4: Relative higher-order effects for different methods: (magenta, long dash dotted) = cr reB , 
(blue/ cyan and dash-dotted/ dashed) = a TCS +, and (red, solid) = <TBom+ vs OBom- 
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